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Abstract 



t-H DAMA annual modulation data and, CoGeNT, CDMS-II, EDELWEISS-II, CRESST excesses of events over the expected back- 
ground are reanalyzed in terms of a dark matter particle signal considering the case of a rotating halo. It is found that DAMA data 
favor the configurations of very high mass dark matter particles in a corotating cold flux. A similar high-mass/low-velocity solution 
^ would be compatible with the observed events in CoGeNT, CDMS-II, EDELWEISS-II and CRESST experiments and could be of 

t^H interest in the light of the positron/electron excess measured by Pamela and Fermi in cosmic rays. 
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1. Introduction 

Since 1996 the sodium iodide experiments of DAMA col- 
laboration (DAMA/Nal and DAM A/LIBRA) have measured an 
annual modulation of the single-hit counting rate which has the 
proper features expected for a dark matter induced signal [ 1 ] . 

More recently, other experiments (CoGeNT fl, CDMS-II 
fiH, EDELWEISS-II 0, CRESST H]) have reported a pre- 
liminary observation of some excess of events relative to the 
expected backgrounds. 

The DAMA annual modulation signal and the other experi- 
ment excesses, if interpreted as dark matter with dominant spin 
independent interaction in the isothermal halo model, implies 
that dark matter particles possess a mass in the range of 5-15 
GeV and an elastic scattering cross section with nucleons in the 
order of 1(T 4 pb HUH. 

In this paper the same data are reanalyzed relaxing the hy- 
pothesis of isothermal halo model, however it is assumed that 
the dark matter local velocity distribution can still be approxi- 
mated as a single Maxwellian flux: 



f(v, v e ) 



3/2' 



(1) 



Here the Earth velocity relative to the dark matter flux is 
given by: v e = v ffi (f) + v - v DM = v ffi (f) + v^ R + v LSR - v~ DM \ 
where: v©(f) is the Earth velocity in the solar system frame; 
jg R = i? - its* = (10.0,5.25,7.17) km/s is the Sun veloc- 
ity relative to the Local Standard of Rest (LSR), and visr - 
(0,220 + 30,0) km/s [10]. Therefore assuming a rotating dark 
matter halo (v DM + 0) one can write: v Ll5S - v DM - (0, Vi ag , 0), 
where V[ ag is the LSR velocity with respect to the dark matter 
flux. Fixing vi ag = vq - 220 km/s the eq. ([1) provides the 
isothermal halo model, however, in this analysis, the vq and vi ag 



parameters are kept free and it is important to note that config- 
urations of vo and vi ag that are far from the isothermal halo ones 
can be physically meaning 

To avoid parameter proliferation, only the case of dominant 
spin independent interaction for elastically scattering dark mat- 
ter will be considered and the effects of uncertainties in the val- 
ues adopted for other parameters (quenching factor, form fac- 
tor, possible presence of channeling, etc..) will be neglected. 
Therefore a four-parameter space (vo,vi ag , Mw,£oo~p) will be 
considered here, where Mw is the particle mass, cr p is the pro- 
ton cross section and £o — icev/cm 3 ' s t ^ le dark matter density 
in units of 0.3GeV/cra 3 . 

2. Experimental observables 

In this section the data used in the analysis are listed for each 
experiment under consideration: 

2.1. DAMA/Nal and DAMA/LIBRA 

The total exposure of 1.17 tonxyr of Nal(Tl) provides three 
complementary observables: 

a) A modulated time behavior in the 2-6 keV window (see data 
in fig. © taken from fig. (4) of ref. |Q]]) 

b) The energy distribution of the observed modulation ampli- 
tude, assuming a fixed phase to — 152.5 d (see data in fig. (HI 
taken from fig. (6) of ref. ll|]). In the following analysis the 
data in the 2-8 keV interval will be considered. 

c) The energy distribution of the unmodulated counting rate 
(see data in fig. (|5]l taken from fig. (27) of ref. 10]). This 
energy distribution provides a limit for the sum of background 
and unmodulated dark matter induced signal and therefore the 
limit of 0.25 cpd/(kg x keV) for the possible unmodulated dark 
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as an example ACDM halo simulations with baryons predict a corotating 
dark disk having vq ~ 50 km/s and vi ag ~ 50 km/s fTltll^l 
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matter induced signal, is cautiously assumed in the following 
analysis; this choice allows large space for the presence of a 
low energy background component in the measured counting 
rate. 

2.2. CoGeNT 

The data of fig. (3) of ref. J2l are considered for the exposure 
of 330g x 56d collected by CoGeNT germanium detector. 

Only the data in the 0.4-0.9 keV window have been used 
for the evaluation of the dark matter allowed configurations in 
the hypothesis that this excess is induced by dark matter elastic 
scattering; however the whole energy interval in fig. (3) of ref. 
J2] is considered for the evaluation of the upper limit. 

2.3. CDMS-II and EDELWEISS-II 

The exposure of 969 kg x d collected by CDMS-II germa- 
nium detectors [4] is considered. Eleven events were observed 
within the recoil acceptance region passing the rejection cuts in 
the 10-150 keV energy range. The neutron background is not 
able to explain the CDMS-II measured events; however some 
of these events could be ascribed to surface background, in par- 
ticular for the low energy region. In the following the hypothe- 
sis that the measured event excess could be due to dark matter 
elastic scattering is considered. Moreover the data of the very 
low energy analysis of CDMS-II (see fig. (1) of ref. 03]) are 
also considered in the evaluation of the upper limit. The re- 
cent result of EDELWEISS-II [5] (where five recoil events are 
measured collecting the exposure of 384 kg x d) seems to be 
compatible with the CDMS-II data; therefore, for simplicity, 
only the CDMS-II data will be considered in this analysis. 

2.4. CRESST 

The preliminary exposure of 564 kg x d collected by 
CRESST - CaWOn detectors is considered [[6j]. In the energy 
range ~ 15-40 keV (the lower threshold is different for dif- 
ferent detectors) 38 events are observed in the Tungsten recoil 
band and 52 events in the Oxygen one. Despite the fact that 
206 Pb recoils from a decay of 210 Po can contribute to the back- 
ground in the higher energy part of the Tungsten recoil band 
and that the Oxygen recoil band is partially overlapped by the 
a recoil band, only a fraction of the observed events can be as- 
cribed to the evaluated background. In this analysis, to account 
for the possible impact of a confirmed excess in CRESST data, 
the case where 30 events of Tungsten recoil and 30 of Oxygen 
recoil are induced by dark matter elastic scattering will be con- 
sidered; this example will be generically addressed as CaW04. 

3. Parameter estimation 

The joint estimation of the four parameters 
(vo, v iag, Mw,%ocr p ) confidence interval has been obtained 
by solving: 

- 2lnL(v , viag, Af w , fto>) + 2lnL max = A (2) 

for the appropriate values of A (A = 7.78 and 13.28 for 90% and 
99% C.L. respectively). In the eq.© L(vo, v/ ag , Mw, £o0>) is the 
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Figure 1: Horizontally hatched areas: allowed configurations (90% and 
99% C.L.) for unconstrained DAMA/Nal + DAMA/LIBRA data. Cross 
hatched areas: allowed configurations (90% and 99% C.L.) for DAMA/Nal 
+ DAMA/LIBRA data combined with CoGeNT, CDMS-II and CRESST data. 
Filled area: configurations having a C.L. better than the one of isothermal halo 
model (t'o = 220 km/s and vi m = 220 km/s) for DAMA/Nal + DAMA/LIBRA 
data unconstrained. 



global likelihood function and L max is the likelihood maximum 
value over the four parameter space. In the estimation of confi- 
dence intervals, the Gaussian approximation has been adopted 
for the likelihood of DAMA/Nal + DAMA/LIBRA data; more- 
over, the possible presence of unknown background explaining 
part or all of the measured events for the other experiments has 
been considered. 

In fig. (Q~|) the projection of the confidence interval surface 
(90% and 99% C.L.) in the plane (vo vs v/ flg ) is shown for two 
cases: 

• unconstrained DAMA/Nal + DAMA/LIBRA data (hori- 
zontally hatched area) 

• DAMA/Nal + DAMA/LIBRA data combined with Co- 
GeNT, CDMS-II and CRESST data (cross-hatched area) 

It can be noted that D AM A/N al + DAM A/LIB R A data favor 
configurations having low velocity dispersion (low Vo) and rel- 
atively low viag which would imply a relatively cold and coro- 
tating dark matter flux in the Galaxy. 

Combining the DAMA data with the constraints from other 
experiments further strengthens this indication. 

As a comparison, in fig. ([TJ, the projection of the config- 
urations having a C.L. better than the one of isothermal halo 
model (vo = 220 km/s and v/ ag = 220 km/s) for unconstrained 
DAMA/Nal + DAMA/LIBRA data, is reported. 
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Figure 2: Allowed configurations in the volume (vo, v; ag , Mw). The confidence 
levels and the cases of unconstrained/constrained DAMA data adopt the same 
palette code used for fig. (TJ. The configurations of very heavy dark matter 
particles are favored. 



In fig. (f2]i the allowed configurations in the volume 
(vo, vi ag , Mw) are shown. The confidence levels and the cases of 
unconstrained/constrained DAMA data adopt the same palette 
code used for fig. It may be noted that configurations hav- 
ing low vo and low V; og require a very high Mw', only a similar 
configuration would provide a flux of dark matter particles with 
enough kinetic energy to allow nuclear recoils events beyond 
the experimental thresholds. 

4. Comparison of the annual modulation signal with re- 
spect to the case of isothermal halo model 

Here, the expected dark matter annual modulation signal fea- 
tures in Nal(Tl) are compared assuming the four different mod- 
els listed in table (Q]). 



Model 


v (km/s) 


v, ag (km/s) 


M w 


f o> (pb) 


a) 


220 


220 


60GeV 


1.3 x lO -5 


b) 


220 


220 


lOGeV 


9.3 x lO" 5 


c) 


10 


95 


90TeV 


5.8 x 10~ 4 


d) 


20 


75 


90TeV 


4.6 x 10~ 4 



Table 1: Models adopted in fig. (5), {4} and (3) 



In fig. (0 the expected modulation behavior in the 2-6 keV 
energy region of DAMA for the four considered models is 
shown. The vertical lines mark the time of \52.5d ~ 2 nd of June 
(where the maximum of the modulation amplitude is expected 
for a non-rotating halo model) and 146 d where a maximum can 
be easily achieved, for example, assuming a corotating flux. 
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Figure 4: Expected energy distribution of modulation amplitudes (S m ) for the 
models listed in table (T}. Data points are taken from fig. (6) of ref. fj]]. 



The data points are taken from fig. (4) of ref. QJJ] and repre- 
sent the annual modulation signal measured by DAMA/Nal + 
DAMA/LIBRA. The measured time of maximum of the modu- 
lation in DAMA is 146+7d, which is compatible both with non- 
rotating as well as with many of the rotating halo models. It is 
important to note that the modulation behavior is roughly sinu- 
soidal but for some extremal models also large departures from 
a pure sinusoid can be found. In fig. © the expected energy 
distribution of modulation amplitudes (S m ) for the four consid- 
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Figure 5: Expected energy distribution of the unmodulated part of the counting 
rate (So) for the four models listed in table {T}. Data points are taken from fig. 
(27)ofref. O. 



Figure 6: Allowed configurations at 2<x C.L. obtained for the isothermal halo 
model with vo = 220 km/s and \>i ag = 220 km/s. 



ered models is shown. The data points are taken from fig. (6) 
of ref. 01 and represent the annual modulation amplitude en- 
ergy distribution measured by DAMA/Nal + DAMA/LIBRA. 
In fig. (0 the expected energy distribution of the unmodulated 
part of the counting rate (So) for the four considered models is 
shown. The data points are taken from fig. (27) of ref. 01311 and 
represent the measured counting rate of DAMA/LIBRA; they 
are the sum of the background and of the possible dark matter 
signal. The dot-dashed line marks the limit of 0.25 cpd/(kg x 
keV) cautiously assumed for the maximum allowed So value in 
this analysis. It is important to note that corotating halo models 
offer a large S m /So ratio allowing the presence of a reasonable 
background component also in the low energy part of DAMA 
data. 

5. Allowed regions fixing the halo: an example. 

In this section, as an example, the halo model will be speci- 
fied to fixed vo and vi ag values. In this fixed halo model, the 2<x 
confidence intervals in the (Mw vs £oo>) plane are evaluated 
considering configurations having A < 6.18 with respect to the 
maximum likelihood of the considered halo model. 

In fig. (0 the allowed regions, assuming the isothermal halo 
model (vo = vi ag = 220 km/s) are shown. 

As a comparison, the allowed regions for the case of a cold 
corotating halo (vo = 20 km/s and v/ flg = 75 km/s) are given in 
fig. 0. In both figures the dashed curve is the limit that can be 
evaluated with CDMS-II when the low energy threshold data 
are also considered [ 14] . It may be noted that compatibility 
among possible positive hints for dark matter could be achieved 
in models of very heavy particles (Mw > few TeV) forming a 
cold corotating halo. 
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Figure 7: Allowed configurations at 2<x C.L. obtained for a cold corotating halo 
with i'o =20 km/s and v/ ag = 75 km/s. The black filled area inside the DAMA 
region marks the configurations allowed at lcr C.L. (A < 2.3). 



As an example, for vo =20 km/s, V[ ag = 75 km/s, Mw = 20 
TeV and £o0> = 10 pb, one would expect: 

• only a fraction of 0.5% of the total CoGeNT rate in the 
0.4-0.9 keV window due to dark matter elastic scattering 

• ~ 13 recoils measured in CDMS-II 

• ~ 60 recoils measured in CaWC<4 (CRESST-like) mainly 
expected to lie in Tungsten band. 
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6. Conclusions 



DAMA annual modulation data and, CoGeNT, CDMS-II, 
EDELWEISS-II, CRESST excesses of events over the expected 
background have been reanalyzed in terms of a dark matter par- 
ticle signal considering the case of a rotating halo. It has been 
found that the data favor the configurations of very high mass 
dark matter particles in a corotating cold flux. This solution is 
intriguing since ACDM halo simulations with baryons predict 
a corotating dark disk having vo ~ 50 km/s and vt ag ~ 50 km/s 

ESQ 

A similar high-mass/low-velocity solution could be of in- 
terest in the light of the positron/electron excess measured by 
Pamela and Fermi in cosmic rays (see e.g. Ill5h . 

Finally, the possibility of multi-component dark matter 
should also be taken into account: one could consider, for ex- 
ample, sterile neutrinos which form a dominant and warmer 
component of dark matter in galaxies together with a sub- 
dominant population of heavy WIMPs mainly gathered in the 
core of dwarf galaxies and, therefore, in the accreted dark disk. 

References 

[1] R. Bemabei et al., Eur. Phys. J. C 67, 39 (2010). 

[2] C.E. Aalseth et al., arXiv: 1002.4703 

[3] Z. Ahmed et al., Science 32 7, 1619 (2010). 

[4] Z. Ahmed et al.. arXiv: 1012.5078 1 

[5] E. Armengaud et al., arXiv: 1 103.4070 

[6] see e.g. Franz Probst slides at the workshop "Dark Matter: Direct Detec- 
tion and Theoretical Developments", Princeton, November 15-16, 2010. 

[7] S. Chang et al., JCAP 1008:018 (2010), arXiv: 1004:0697. 

[8] A. Bottino et al., P hys. Rev. D 81, 10 7302 (2010). 

[9] N. Fornengo et a UarXiv: 101 1.4743 1 

[10] C. McCabe, Phys. Rev. D 8 2, 023530 (2010). 

[11] J. Read et al.. larXiv:0803.27T4l 

[12] C.W. Purcell et al., Astrophys J. 703 2275 (2009). 

[13] R. Bernabei et a l., Nucl. Ins tr. and Meth. A 592, 297 (2008). 

[14] Z. Ahmed et al.. larXiv:1011.2482l 

[15] I. Cholis and L. Googenhough, JCAP1009:010 (2010), arXiv: 1006.2089 



5 



